Future constraints on the Hu-Sawicki modified gravity scenario 



Matteo Martinelli^ , Alessandro Melchiorri^, Olga Mena^, Valentina Salvatelli^ and Zahara Girones^. 
^Physics Department and INFN, Universita' di Roma "La Sapienza", Pie Aldo Mora 2, 00185, Rome, Italy, and 
^ IFIC-CSIC and Umversidad de Valencia, Valencia, Spam 
(Dated: September 23, 2011) 

We present current and future constraints on the Hu and Sawicki modified gravity scenario. 
This model can reproduce a late time accelerated universe and evade solar system constraints. 
While current cosmological data still allows for distinctive deviations from the cosmological constant 
picture, future measurements of the growth of structure combined with Supernova la luminosity 
distance data will greatly improve present constraints. 
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I. INTRODUCTION 



Astrophysical observations have led to the inference 
that our universe is approximately flat and its mass- 
energy budget consists of 5% ordinary matter, 22% non- 
baryonic dark matter, plus a dominant negative-pressure 
component that accelerates the Hubble expansion 
The current accelerated expansion can be explained as 
the presence of a Cosmological Constant (CC hereafter) 
associated to the energy of the vacuum in Einstein's equa- 
tions. However, the naively theoretically expected value 
exceeds the measured one by 123 orders of magnitude 
and it needs to be cancelled by extreme fine-tuning. A 
dynamical alternative attributes the accelerated expan- 
sion to a cosmic scalar field, quintessence (5l-[To| in which 
case the equation of state w could vary over time. How- 
ever, these quintessence models are not better than the 
CC scenario as regards fine-tuning. 

In this paper we focus on the third possible scenario, 
in which the gravitational sector is modified, as an alter- 
native to exp lain the observed cosmic acceleration (see 
Refs. |ll| - [l3| and references therein). Although this re- 
quires the modification of Einstein's equations of gravity 
on very large distances [l4|, or on small curvatures |l5l - 
[iTj . this is not unexpected for an effective 4-dimensional 
description of higher dimensional theories. Modifications 
of gravity have been examined in the context of accel- 
erated expansion. The proposed modified gravity mod- 
els have extra spatial dimensions or an action which is 
non linear in the curvature scalar, that is, these mod- 
els include extensions of the Einstein-Hilbert action, for 
instance, to higher derivative theories 18], scalar-tensor 
theories or generalized functions of the Ricci scalar f{R). 
Modified gravity models have been confronted with cur- 
rent and future data extensively in the literature [Tol - 
Issj . However, it is well known that f{R) gravity models 
that produced late time acceleration also have problems 
to pass solar system tests [36l - l46j . The reason is that 
f{R) gravity theories introduce a scalar degree of freedom 
given by fn. = df /dR that, for the background cosmolog- 
ical density, is very light. As a consequence, it produces 
a long-range fifth force, leading to a dissociation of the 
space-time curvature from the local density. Then, the 
metric around the sun is predicted to be different than 
what is observed. Chiba [43 has shown the conditions 



under which a given /(i?) model is equivalent to a scalar- 
tensor theory with Parametrized Post-Newtonian (PPN) 
parameter 7 = 1/2, far outside the range allowed by ob- 
servations, I7 — 1| < 2.3 • 10~^ [15. However, some f{R) 
theories are still viable: the scalar field mass could be 
large and therefore it would not have an effect at solar 
system scales. Another possibility is a scale dependent 
scalar field mass, as in the chameleon mechanism [49- 53,]. 
In chameleon cosmologies, the effective mass of the scalar 
field becomes very large in high density environments (as 
in the Sun's interior) and the induced fifth force range 
would be below the detectability level of gravitational 
experiments. 

Among a plethora of f{R) models, we focus here on the 
one proposed by Hu and Sawicki [42*1 (HS hereafter). This 
model is able to reproduce the late time accelerated uni- 
verse, but with distinctive deviations from a cosmological 
constant. The model has also been shown to satisfy the 
conditions needed to produce a cosmologically viable ex- 
pansion [s^l- More interestingly, the model is designed 
to posses a chameleon mechanism that allows to easily 
evade solar system constraints. The authors of Ref. |54| 
have analyzed the HS model exploiting current cosmolog- 
ical data. Here we update these results using the most 
recent cosmological data and we present the expected 
constraints on the HS model from future Baryon Acous- 
tic Oscillation (BAO) surveys and future Supernovae la 
(SNIa) luminosity distance measurements. 

The structure of the paper is as follows. In Sec. HIl we 
briefly describe the HS model. In Sec. lIIII we present con- 
straints on this model using the most recent cosmological 
data. Section ITVl describes the method used here to fore- 
cast data from future BAO and SNIa surveys. Future 
constraints from the former data are presented in SecfVl 
We draw our conclusions in Sec. IVIl 



II. THE MODEL 

We briefly review below the basic equations and results 
of the HS model. This model has a modifled Einstein- 
Hilbert action: 
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where Cm is the matter lagrangian, = SttG and 



(2) 



with = K^po/S, being po the average density today 
and ci, C2 and n as free parameters. 

Varying the action Eq. ([T]) with respect to the metric 
g'^'^ one obtains the modified Einstein equations 

G^. + - - afj,)g,,, - V^V.fR = k^T^. , (3) 

where fa — df /dR and /rr = d^ f /dE? . Assuming a fiat 
Friedmann Robertson Walker (FRW) metric, the modi- 
fied Friedmann equation reads 



_ f^i^HH' + H^) + {+ H^JbrR' = ^ , (4) 
6 3 



with ' = d/dlna. Defining the new variables yn — 
{H'^ /m?) — and yji — {R/m^) — 3a~^ the Friedmann 
equation can be written as a system of two ordinary dif- 
ferential equations: 
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Using these last two equations we can relate ci and 
C2 to the free parameters of the model, n and fu^ = 
fnilna = 0), and to Qm- The parameter is con- 
strained to l/flol < 0.1 by solar system tests [i^l and 
therefore we will not investigate larger values in the next 
sections. Figures [T] and [2] show the behavior of w as a 
function of the redshift z for different values of /^jg and 
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In order to compare the HS model with the cosmolog- 
ical constraints usually derived under the assumption of 
a dark energy fluid, it is useful to introduce an effective 
dark energy component with a present energy density 
rix = 1 — rirn and a equation of state w, where flm is 
the effective matter energy density at present time. Of 
course, in reality, no dark energy component is present 
and the only component of the universe is matter, be- 
ing modified gravity the responsible for the accelerated 
expansion. Considering the Friedmann equation 
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the effective equation of state parameter w for the dark 
energy component is given by 
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The free parameters ci and C2 that appear in Eq. ^ 
can be expressed as a function of the effective density 
parameters 



FIG. 1. The equation of state parameter w for n = 1 and 
/flf, — —0.1, —0.03 and —0.01 (solid, dashed and dotted lines). 
Q,m is set to 0.24. 

As we can see in FigH] and Figl2] the equation of state 
parameter w follows a peculiar behavior as a function 
of the redshift. At the present time (z = 0) w has al- 
ways a value higher than the one predicted by the ACDM 
model {w — —I) and, moving towards higher redshifts, 
it decreases towards the phantom region, i.e., taking val- 
ues lower than —1. For even higher redshifts, w moves 
asymptotically towards —1. If the absolute value of fa^ 
is decreased, w gets closer to —1, while if the parame- 
ter n is increased the phantom crossing occurs at lower 
redshifts. 

After solving the background equations in terms of 
the new variables yn = [H'^/m?) — a^^ and yn = 
[R/m?) — 2>a~^ we are able to predict the expected theo- 
retical values for a set of observables in order to constrain 
the free parameters of the HS model [n and /aq)- We 
consider here the luminosity distance 
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FIG. 2. The equation of state parameter w for n = 2 and 
/rq equal to —0.1 (solid line), —0.03 (dashed line) and —0.01 
(dotted line). Q,m is set to 0.24. 



as well as the Hubble parameter 



distances. For the BAO, we exploit data from the Sloan 
Digital Sky Survey ^ and from 2-Degree Field (2dF) 
Galaxy Redshift Survey In addition, data from the 
WMAPril] satellite has been used to extract the CMB 
derived parameters R (the shift parameter) and I a (pa- 
rameter related to the first peak position). We build a 
global variable 

(16) 
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The likelihood function is defined as 
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and is marginalized over the present value for the Hubble 
expansion rate Hq. The free parameters of the HS model 
are n and fuQ. The value of the present dark matter en- 
ergy density flm has been fixed to 0.27. Figure |3] shows 
the constraints on the HS model. Notice that current 
background data prefer small values of both n and fm- 
These bounds are stronger than those presented in [54| . 
where R and Ia were not used and the SNIa and BAO 
datasets were from older catalogues. In the following sec- 



H{a) = ^n„,H-i{yH + a-3) , 
and the angular diameter distance 



dA^a) = 



^ da 1 /"^ da 



a a'^H{a) Ho 7„ 



(12) 



a'^y^miVH + a ^) 



(13) 

We also compute here the growth of structure pre- 
dicted by the HS model. The linear growth equa- 
tion for modified gravity models is scale dependent and 
reads 0, HI 
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where ' = d/da, d is normalized such that 6 ^ a when 
a — > and the factor Q is given b 
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Note that in general relativity Q is zero and therefore 
the linear density growth is scale independent for all dark 
energy models. However, for f{R) models, the scale de- 
pendent Q{k, a) induces a nontrivial scale dependence of 
the growth 5. 



FIG. 3. Constraints on n and /_ro using current data. 

tions we show that future growth data will prefer higher 
value of n values and the tension between future growth 
and background data will tighten the constraints on the 
HS model. 



III. CURRENT CONSTRAINTS 

We exploit here current data from standard candles 
(SNIa) and standard rulers (BAO, CMB) in order to con- 
strain the expansion history in the HS modified gravity 
model. Our analysis make use of the Union2 survey [j] , 
which provides 557 SNIa useful to constrain luminosity 



IV. FUTURE DATA 

A. SN Luminosity distance data 

We shall exploit future SNIa data on the reduced mag- 
nitude 
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being the luminosity distance given by Eq. (|TT|) . We 
consider here a mock catalog of 2,298 SNIa, with 300 
SNIa uniformly distributed out to z — 0.1, as expected 
from ground-based low redshift samples, and an addi- 
tional 1998 SNIa binned in 32 redshift bins in the range 
0.1 < z < 1/7, as expected from JDEM or similar fu- 
ture surveys |57|. We have considered both intrinsic and 
systematic errors on the reduced magnitude. 
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FIG. 4. Forecasted SNIa luminosity distance data for a 
JDEM-like survey. 



is the cosine of the angle between the line of sight and 
the wavevector fc. Notice that perturbations with fc per- 
pendicular to the line of sight are not distorted. The 
relation among real space and redshift space overdensi- 
ties given by Eq. (|T9l) was first derived by Kaiser 
and it arises from the continuity equation, which relates 
the divergence of the peculiar velocity to the dark matter 
overdensity 8. Redshift space distortions, then, relate pe- 
culiar velocities to the logarithmic derivative of the linear 
growth factor, /. A measurement of /3 = //& will provide 
information on the growth of structure formation if the 
galaxy bias h is known. 

In our analysis we have computed the forecasted errors 
for / assuming a ACDM fiducial model by means of a 
Fisher matrix analysis, marginalizing over the bias h. We 
focus here on the BAO experiments BOSS [59] and Euclid 
[60j . For the BOSS (Euclid) experiment we assume six 
(nineteen) redshift bins ranging from z — 0.15 to z = 0.65 
{z = 0.15 to z = 1.95) and a galaxy survey area of 10000 
(20000) deg^. The mean galaxy densities for these two 
experiments are considered to be constant with values of 
2.66 X lO"'' and 1.56 x lO'^ h Mpc'^ for the BOSS and 
the Euclid surveys respectively. 

We have combined t he g alaxy survey fisher matrices 



with the CMB Planck 61 



Fisher matrix, see Ref [62| . 
The expected errors on / are depicted in Figs. [5] and [S] 
and are similar to the ones obtained using the numerical 
tools of Ref. [Hi. 



In the following, to produce mock SNIa luminosity dis- 
tance and mock growth data from BAO surveys we shall 
assume a ACDM model with the following cosological pa- 
rameters: rjfc/i^ = 0.02258, ^Ji^ = 0.1109, = 0.963, 
T = 0.088, As = 2.43 x 10^^ and 9 ^ 1.0388, which 
correspond to the best-fit values from the WMAP seven 
year data analysis, see Ref. |l|. 



B. Growth data 

Galaxy surveys measure the redshift of the galaxies, 
providing, therefore, the redshift space galaxy distribu- 
tions. From those redshifts the radial position of the 
galaxies are extracted. However, the inferred galaxy dis- 
tribution (and, consequently, the power spectrum) is dis- 
torted with respect to the true galaxy distribution, be- 
cause in redshift space one neglects the peculiar velocities 
of the galaxies. These are the so called redshift space dis- 
tortions. 

In linear theory and with a local linear galaxy bias b 
the relation between the true spectrum in real space and 
the spectrum in redshift space reads 
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where /3 = f/b, being / the logarithmic derivative of 
the linear growth factor S{a) given by Eq. and /ifc 




0,4 0,5 0,6 0,7 



FIG. 5. Forecasted BOSS data for /, the logarithmic deriva- 
tive of the linear growth factor, as a function of the redshift. 



V. FUTURE DATA ANALYSIS AND RESULTS 

In this section we perform a numerical fit to the HS 
model using the mock data sets described above by means 
of a analysis. We compute two different functions. 
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FIG. 6. Forecasted Euclid data for /, the logarithmic deriva- 
tive of the linear growth factor, as a function of the redshift. 



one associated to growth of structure (xlrowth) 
other related to SNIa luminosity distance data, Xsn- 
final likelihood will exclusively depend on the free param- 
eter of the HS model, n and ffjg , thus for the two different 
analyses (xlrowth ^^"^ Xsn) ^^'^ ^'^^^ f*-"^ ^^"^ com- 



bined one ixtot = X growth + Xsn) we marginalize over 
flm- For the growth of structure we compute the xirowth 
at different scales k scales and perform a mean over k 
from k = O.Ol/i Mpc'^ to fc = O.lh Mpc-\ 
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FIG. 7. 1 and 2-cr contours in the (n, /hq) plane arising from 
SNIa measurements of the luminosity distance (black solid 
and dashed lines) and those arising from the BOSS measure- 
ments of the linear growth of structure (filled regions). 



FIG. 8. 1 and 2-cr contours in the (n, fn^^) plane arising from 
SNIa measurements of the luminosity distance (black solid 
and dashed lines) and those arising from the Euclid measure- 
ments of the linear growth of structure (filled regions). 



Figure [7] shows the 1 and 2-cr contours on the (n, fji^) 
plane arising from the analysis of SNIa luminosity dis- 
tance data as well as from measurements of the linear 
growth of structure from the BOSS experiment. The 
empty regions depicted by the dotted and solid lines show 
the results from SNIa luminosity distance data, while the 
filled regions depict the results from the measurements of 
the linear growth of structure. Notice that background 
and growth measurements are complementary, allowing 
different regions in the parameter space. The combina- 
tion of both data sets provide a very powerful tool to 
severely constrain the HS model. 

Figure [8] shows the equivalent to Fig.[7]but considering 
measurements of the linear growth of structure from the 
Euclid galaxy survey. Notice that the currently allowed 
region will be ruled out if the true cosmology is a ACDM 
scenario. Therefore, the BOSS and Euclid surveys, com- 
bined with background SNIa luminosity distance mea- 
surements, will be able to exclude a very large region of 
the parameter space of the HS model, allowing only for 
very high values of n for which ci/c^ — ?► 0, see Figs. [9] 
and [10] in which the combined analyses xj^f are shown. 
Figure fTTl illustrates that, in the limit of ci/c^ — > 0, i.e. 
for very large values of n, f{R) is equivalent to a cosmo- 
logical constant. Therefore, future data from the Euclid 
galaxy survey combined with SNIa measurements will be 
able to recover the true fiducial cosmology even if the 
data is fitted to a non constant f(R). 




n 



FIG. 9. 1 and 2-a contours in the (n, /hq) plane arising from 
BOSS measurements of the hnear growth of structure com- 
bined with SNIa luminosity distance data from the JDEM 
survey. Note that the axis ranges are different to those of 
Fig. El 
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FIG. 10. 1 and 2-<j contours in the {n, fng) plane arising 
from Euclid measurements of the linear growth of structure 
combined with SNIa luminosity distance data from the JDEM 
survey. Note that the axis ranges are different to those of 
Figs. □] and H 
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FIG. 11. The top and bottom panels depict the effective equa- 
tion of state w{z) for n = 9.8 and fug = —0.012, values that 
are in agreement with Euclid growth data. The top panel 
shows w{z) using the same ranges than those use in Figs. [1] 
and [2] The bottom panel shows an axes range that allows to 
appreciate the behaviour of w{z). 



data allow for a region of the parameter space where the 
parameter n of the model can have small values, induc- 
ing a late time dynamical behavior on the effective dark 
energy equation of state w{z). Future constraints on the 
Hu-Sawicki modified gravity model arising from measure- 
ments of the linear growth of structure and of SNIa lu- 
minosity distances are also presented. While luminosity 
distance data allow for small values of n, the growth of 
structure data prefer higher values of n. The combina- 
tion of these two observables allows to tightly constrain 
the HS model. We have generated mock growth and lu- 
minosity distance data data for a fiducial ACDM model 
and fitted these data in a Hu-Sawicki modified gravity 
scenario. The reconstructed effective dark energy equa- 
tion of state almost is identical to that of a ACDM model. 



VI. CONCLUSIONS 



The Hu-Sawicki modified gravity scenario is analyzed 
with current background expansion data, updating pre- 
vious results [53|. Our analysis shows that background 
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